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' Abstract 

' _ j ' It is shown that a spacetime with coUisionless matter evolving from 
' data on a compact Cauchy surface with hyperbohc symmetry can be glob- 
es ■ ally covered by compact hypersurfaces on which the mean curvature is 

constant and by compact hypersurfaces on which the area radius is con- 

I I stant. Results for the related cases of spherical and plane symmetry are 

bX). reviewed and extended. The prospects of using the global time coordinates 

^ ' obtained in this way to investigate the global geometry of the spacetimes 

. 1-H , concerned are discussed. 

X' 

H 

5r ! 1 Introduction 

One of the key mathematical problems in general relativity is the determination 
of the global properties of solutions of the Einstein equations coupled to various 
matter fields. We investigate this problem in the case of coUisionless matter 
described by the Vlasov equation. The underlying strategy is to first estab- 
lish the existence of a suitable global time coordinate t and then to study the 
asymptotic behaviour of the solution when t tends to its limiting values, which 
might correspond to the approach to a singularity or to a phase of unending 
expansion. The following is almost exclusively concerned with the first of these 
two steps. Since the general case is beyond the range of current mathematical 
techniques we study a case with symmetry. 

There are several existing results on global time coordinates for solutions of 
the Einstein- Vlasov system. In the spatially homogeneous case it is natural to 
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choose a Gaussian time coordinate based on a homogeneous hypersurface. The 
maximal range of a Gaussian time coordinate in a globally hyperbolic solution 
of the Einstein- Vlasov system evolving from data on a compact manifold which 
are homogeneous (or locally homogeneous) was determined in [ pO) . It is finite 
for models of Bianchi IX and Kantowski-Sachs types and finite in one time di- 
rection and infinite in the other for the other Bianchi types. All other results 
presently available on the subject concern spacetimes which admit a group of 
isometries acting on two-dimensional spacelike orbits, at least after passing to a 
covering manifold. The group may be two-dimensional (local U{1) x U{1) sym- 
metry) or three-dimensional (spherical, plane or hyperbolic symmetry). In all 
these cases the quotient of spacetime by the symmetry group has the structure 
of a two-dimensional Lorentzian manifold Q. The orbits of the group action 
(or appropriate quotients in the case of a local symmetry) are called surfaces 
of symmetry. Thus there is a one-to-one correspondence between surfaces of 
symmetry and points of Q. 

Three types of time coordinates which have been studied in the inhomoge- 
neous case are CMC, areal and conformal coordinates. A CMC time coordinate 
t is one where each hypersurface of constant time has constant mean curvature 
and on each hypersurface of this kind the value of t is the mean curvature of 
that slice. In the case of areal coordinates the time coordinate is a function of 
the area of the surfaces of symmetry. In some papers in the literature it is taken 
to be proportional to the area. In this paper it is taken to be proportional to 
the square root of the area. In the case of conformal coordinates the metric on 
the quotient manifold Q is conformally flat. The properties of the last two kinds 
of time coordinates will be described in more detail later. 

Next the known results concerning these time coordinates in solutions of the 
Einstein- Vlasov system will be summarized briefly. For the detailed statements 
the reader is referred to the original papers. In the case of spherical symmetry 
the existence of one compact CMC hypersurface implies that the whole space- 
time can be covered by a CMC time coordinate which takes all real values ||l9|, ^ . 
The existence of one compact CMC hypersurface in this case was proved later 
by Henkel , thus providing a complete picture in the spherically symmetric 
case. His proofs use the concept of prescribed mean curvature (PMC) folia- 
tion which may be useful more generally. A general local existence theorem 
for PMC foliations was proved in In the case of U{1) x £7(1) symmetry 

(which includes plane symmetry as a special case) it was proved in that the 
existence of one compact CMC hypersurface (without loss of generality with 
negative mean curvature) implies that there is a foliation by compact CMC 
hypersurfaces where the mean curvature takes on all negative values. It was 
also shown that the foliation covers the whole spacetime between the initial hy- 
persurface and the singularity, but it was left open whether it covers the whole 
spacetime in the other time direction. The latter question will be discussed fur- 
ther in Section |^. In the special case of Gowdy symmetry (which includes plane 
symmetry) the existence of one compact CMC hypersurface follows from . 
Finally, in the case of hyperbolic symmetry, the results of jl9| imply that the 
existence of one compact CMC hypersurface leads to the same conclusions as in 
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the U{1) X ?7(l)-syinnietric case, provided the Hawking mass is positive on the 
initial hypersurface. (The question of whether the CMC fohation includes all 
negative numbers as values of the mean curvature was not discussed in [|l9j ; it 
will be settled in Section ^.) Under the same positivity condition the existence 
of one compact CMC hypersurface is proved in . 

The main concern of this paper is with the areal and conformal time coordi- 
nates. In fact it is the first of these which is of fundamental interest — conformal 
coordinates serve as a convenient tool in an intermediate step in the proofs. 
Note that the areal time coordinate is automatically positive so that the largest 
possible time interval on which the solution could exist is ]0,oo[. It was shown 
in [p^ for spacetimes with spherical symmetry that if there is one symmetric 
Cauchy surface of constant areal time and if the data on that hypersurface sat- 
isfy a certain inequality then the past of the initial hypersurface is covered by 
an areal time coordinate and this coordinate takes on all values in the range 
]0,i?i], where i?i is its value on the initial hypersurface. A result improving this 
time interval cannot be expected since these spacetimes tend to recoUapse, and 
then the areal coordinate will break down at some point in the time direction 
corresponding initially to expansion. The results on the areal coordinate in |^ 
for the case of plane symmetry are superseded by those on Gowdy symmetry in 

@- 

For Gowdy symmetry the spacetime is globally covered by an areal time 
coordinate which takes all values in the range ]i?o,oo[ for some Ro>0- In some 
cases, such as the vacuum Gowdy spacetimes, it is known that i?o = 0, cf. Q, 
but it is not understood in what generality this holds. For hyperbolic symmetry 
a direct analogue of the result for restricted data in spherical symmetry stated 
above was also proved in . 

In the case of hyperbolic symmetry two main questions were left open by 
the work up to now. Firstly, is there always one symmetric Cauchy surface of 
constant areal time? Secondly, is there an areal time coordinate whose level 
surfaces are compact and symmetric which covers the whole spacetime, and 
does it take on arbitrarily large values? These questions will be answered in 
the affirmative in the following (Theorem 6T). The proofs are modelled on the 
approach of which in turn was inspired by the work of 1^ on the vacuum 
case. A conformal time coordinate is used in an intermediate step. 

In Section 2 the definitions of hyperbolic symmetry and the Einstein- Vlasov 
system are recalled. The local existence of conformal coordinates in a neigh- 
bourhood of the initial hypersurface is demonstrated. In Section 3 the evolution 
of the spacetime in the expanding direction is analysed. A global existence theo- 
rem in areal time is proved which answers the second main question modulo the 
first one. In the section after that a long-time existence result in the contracting 
direction is proved using conformal coordinates. In Section 5 the existence of 
a symmetric compact Cauchy surface of constant area radius is demonstrated, 
thus answering the first main question. All these results are combined to produce 
the main results of the paper in Section 6. Theorem 6.1 asserts the existence 
of a global foliation by hypersurfaces of constant areal time in solutions of the 
Einstein- Vlasov system with hyperbolic symmetry, and Theorem 3.2 makes an 
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analogous statement for hypersurfaces of constant mean curvature. The last 
section discusses possible extensions of these results and looks at what is known 
about the asymptotic behaviour of solutions. 

2 Hyperbolic symmetry 

In figl a definition of spacetimes with surface symmetry was given. This com- 
prised three cases, namely spherical, plane and hyperbolic symmetry. In the fol- 
lowing the case of most interest is that of hyperbolic symmetry. The spacetime 
M is diffeomorphic to IR x S*^ x where F is a compact orientable surface of 
genus greater than one. The manifold M has a covering space M diffeomorphic 
to IR X S'^ X _F with projection p induced by the projection p from the universal 
cover F to F according to p{x,y,z) = {x,y,p{z)). The surface F admits metrics 
of constant negative curvature, and F endowed with the pull-back of any one of 
these metrics is isometric to the hyperbolic plane. Let G be the identity com- 
ponent of the isometry group of the hyperbolic plane. If F is identified with the 
hyperbolic plane then an action (j) oi G on F \s obtained. Define an associated 
action on M by 4>(x,y,z) = {x,y,4>{z)). A spacetime with underlying manifold 
M defined by a metric gap and matter fields is said to have hyperbolic symme- 
try if the pull-back gap of the metric and the pull-back of the matter fields to 
M via p are invariant under (p. The surfaces in M diffeomorphic to F defined 
by the product decomposition will be called surfaces of symmetry. A Cauchy 
surface S will be called symmetric if it is a union of surfaces of symmetry. Then 
S = p~^{S) is invariant under the action of the group G. It is now clear how to 
define abstract Cauchy data for the Einstein-matter equations with hyperbolic 
symmetry. They should be defined on a manifold S of the form x F and 
should be such that the pull-back of the data to S — S^xF under the natural 
covering map is invariant under the natural action of G. The quotient of S by 
the action of G is diffeomorphic to S^. 

Consider now a choice of matter fields for which the Cauchy problem for the 
Einstein-matter equations is well-posed. The example of interest in the following 
is that of collisionless matter satisfying the Vlasov equation. Corresponding to 
initial data for the Einstein-matter equations with hyperbolic symmetry there 
are data on the covering manifold S which have a maximal Cauchy develop- 
ment on a manifold M. We will now construct a certain coordinate system on 
a neighbourhood of the initial hypersurface in M. As remarked in p!9| coordi- 
nates can be introduced on the initial hypersurface S so that the metric takes 
the form A'^{9)d9'^ + B'^{9)d'E'^ , where dT,'^ — dx'^ + sinh^ x dy"^ . Correspondingly 
there are local coordinates on S (local in x and y) where the metric takes this 
form. On a neighbourhood of the initial hypersurface in M we can introduce 
Gauss coordinates based on the given coordinates on S. The group G acts as 
a symmetry group on the initial data. On general grounds the action of G on 
S extends uniquely to an action on the maximal Cauchy development M by 
symmetries (see e.g. Q). The hypersurfaces of constant Gaussian time are in- 
variant under the action of G, as are the hypersurfaces of constant 9. Using the 
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isotropy group of any point it can be seen that the 9x and 9y components of the 
metric of these hypersurfaces must vanish. In the same way, the restriction of 
the second fundamental form of the hypersurfaces to the surfaces of symmetry 
must be proportional to the metric with a factor depending only on t and 9. By 
integration in t it then follows that the spacetime metric takes the form 

(is2 = _dt^+ {t, e)de'' + (t, e)dY? (2. i) 

in the region covered by Gauss coordinates. From this we see that we can 
form the quotient of this metric by the action of a discrete group to get a 
spacetime on a subset of M which is a Cauchy development of the original data 
on S. Furthermore we can form the quotient by G to get the Lorentz metric 
—dt^ + A^{t,9)dO'^ on a subset Q of IRx 5*^ referred to in the introduction. On 
Q we can pass to double null coordinates {u,v) on a neighbourhood of the 
quotient of the initial hypcrsurface. Defining new coordinates by t = ^{u — v) 
and 6 = ^{u + v) puts the metric on Q in conformally flat form. By pull-back 
these define new coordinates on M where the metric takes the form 

ds2 ^ e^r,i^_dt^ + de^) + R^idx"^ + sinh^ xdy^) (2.2) 

where rj and R are functions of t and 9 which are periodic in 9. This proves the 
existence of a conformal coordinate system close to the initial hypersurface. It 
is possible to choose the double null coordinates in such a way that the initial 
hypersurface coincides with t = and the period of the functions rj and R is one. 

To conclude this section we formulate the Einstein- Vlasov system which 
governs the time evolution of a self-gravitating coUisionlcss gas in the context 
of general relativity; for the moment we do not assume any symmetry of the 
spacetime. All the particles in the gas are assumed to have the same rest mass, 
normalized to unity, and to move forward in time so that their number density 
/ is a non-negative function supported on the mass shell 

PM:={g^0p"p^ = -l,p°>O}, 

a submanifold of the tangent bundle TM of the space-time manifold M with 
metric gaf)- We use coordinates {t,x°') with zero shift and corresponding canon- 
ical momenta p"; Greek indices always run from to 3, and Latin ones from 1 
to 3. On the mass shell PM the variable becomes a function of the remaining 
variables {t,x°-,p''): 

The Einstein- Vlasov system now reads 

J -Po 
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where F^^ are the ChristofFel symbols, \g\ denotes the determinant of the metric, 
G"^ the Einstein tensor, and T"'^ is the energy-momentum tensor. The puU- 
back of the number density / to TM is assumed to be invariant under the action 
induced on TM by the action of G on M, a fact which can be used to reduce 
the number of independent variables. 

3 The expanding direction 

In this section we want to investigate the Einstein- Vlasov system with hyperbolic 
symmetry in the expanding direction. We write the system in areal coordinates, 
i.e., the coordinates are chosen such that R = t. The circumstances under which 
coordinates of this type exist are discussed later. We prove that for initial 
data on a hypersurface of constant time corresponding solutions exist for all 
future time with respect to the areal time coordinate. It should be noted that 
our time coordinate has the geometric meaning of the curvature radius of the 
hyperbolic spaces which form the orbits of the symmetry action. Since it requires 
little additional effort we include for the sake of comparison the case of plane 
symmetry and write the metric in the form 

ds2 ^ „e2^dt2 + e^^d9^ + [dx^ + sixvlxdy"^), (3.1) 

where /x and A are functions of t and ^, periodic in 9 with period 1, and 

{sina; for e = l, 
1 for e = 0, 
sinhcc for e~—l. 

For the case e = 1 the orbits of the symmetry action are two-dimensional spheres. 
In this spherically symmetric case the global result below is easily seen to be 
false, cf. so this case will not be considered further. For the case e = the 
orbits of the symmetry action are flat tori, and the coordinates x and y range 
in the interval [0,27r]. 

Lemma 3.1 Assume that (the pull-hack of) f is invariant under the group 
action on TM associated to hyperbolic symmetry. Then f depends only on 

t, 9, p\ (p2)2 + sinh2x(/)2. 

Actually, we will write 

f^f{t,9,w,L) 

where 

w:=e^p\ L -.^ t\ip^ f + sinh^ xip^f). 

In these variables 

p°^e-^'y/^ + w^ + L/t^=■.e-^'(p), 

and L is a conserved quantity along particle orbits. The analogous result holds 
in the case of plane symmetry e = 0. 
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Proof of Lemma 3.1. To see the above assertion on the form of / we proceed as 
follows. In the coordinates 

— =tcoshx, a;^ = tsinhscos?/, a;'^ = tsinhxsin?/ 

on the covering manifold AI the symmetry orbits are given as {x'^ = const} x 
H'^{t) where 

H^t) := {{x\x^,x^) e IR3 I > 0, ~{x^f + {x^ + (5^)2 = _t2}^ 

a hyperbolic space of radius t, cf. p. 108 ff]. The reason for using these 
coordinates is that the isometry group of the hyperbolic spaces now has a nice 
representation: Its elements are given by the restriction to H'^{t) of the linear 
maps X — [x^ ^xP' ^x^Y'' ^ Sx with matrices S € GL(3,IR) which have the property 
that 

{Su,Sv) — {u,v) , u,v£jR^ 

where 

{u,v) ■.= -u'^v'^ +u'^v'^ + u^v-^ 

is the usual inner product on three-dimensional Minkowski space. In addition, 
the matrices S are required to preserve orientation, i. e., 

g^fax b \ ay>o,a5eIR2x2 ^^^^ detas>0, 5,ceIR^ 
\ c as J 

That these matrices represent the isometrics of the hyperbolic spaces can be 
found in ||l^, p. 239 ff] or can be seen directly: In the new coordinates 

- 2A - n - , /-I 2!i\ ^ac^ ^hd^ 
500=e , 50a = 0, 3ab = ea6+(l-e ^) — 7a6 

where 

eil = -l, €22=633 = 1, eab = 0, a ^ 6, 

and 

7ab = 76a = -l, a = l,& = 2,3, 7afc = l othcrwisc, 

and since for the matrices introduced above one has —SiaSih + S2aS2b + 
S'iaS'ib = ^ab, i- e., their columns are orthonormal with respect to (•,•), their 
isometry property is easily seen. Now denote by the canonical momenta 
corresponding to the new coordinates a;", i. e., 

in particular fp=p^. Fix some 9 = xP ,t>Q, x€H^{t), p°GlR, and p = 
(p\p2^p3yr g j^3^ -^g^ g isometry such that 5'x = (t,0,0)*'' and T a ro- 

tation in the [x^ ,x^)-pla.ne such that TSp^ {q^ ,q'^ ,0)^" with q'^>0. Since TS 
represents an isometry under which / should be invariant, 

fi9,x,p°,p) = f{9,TSx,p\TSp)^f{9,t,0Ap\q\q^0) 
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and ^ ^ 

q^^--{TSx,TSp)^--{x,p), 



3\2 



When these quantities are re-expressed in the old momentum variables via (3.2) 
we find 

which proves the assertion in the hyperbolic case. For the plane case the argu- 
ment is even easier. 

In the variables which we have now introduced the complete Einstein- Vlasov 
system reads as follows: 

^tf + ^-:^^ef~{\tw + e^^-^^xg{p))^^f = Q, (3.3) 



e-2^(2Ut + l) + e = 87rtV, (3-4) 

e-^^'{2tfit - 1) - e = STTi^p, (3.5) 

|I0 = -i^:te^'+^j, (3.6) 

(Mee+Me(/^e-Ae)) -e-2A'(^Att-K(At + l/t)(At-/it)) ^Anq, (3.7) 



where 



/CO f>00 
/ {p)f{t,e,w,L)dLdw^e-^''To„{t,0), (3.8) 
-oo Jo 

p{t,e):^-j —fit,9,w,L)dLdw = e-^^T^,it,9), (3.9) 



j{t,0):=^l I wf{t,9,w,L)dLdw = ~e^+''Toi{t,e), (3.10) 



-oo JO 

fOO POO T o 

TT/ / L , , „ ,,,,, I 



q{t,9):^^J —f{t,e,w,L)dLdw^-^T22{t,e)- (3.11) 



recall that 



{p) = v/l+U-2+L/f2. 

We prescribe initial data at some time i = to > 0, 

fito,9,w,L)=fi9,w,L), X{to,9) = \{e), ti{to,9)=fi{9) 

and want to show that the corresponding solution exists for all t>to. To this 
end we make use of the continuation criterion in the following local existence 
result: 
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Proposition 3.2 Lei f €C^(JR'^ x]R^) with f{e + l,w,L)=f{9,w,L) fo 
(6l,w,L)e]R^x]R+, f>0, and 

sup I \w\ I {0,w,L) e supp/| < oo, 
sup |l I {6,w,L) e supp/| < oo. 

Let XjeC^iWi) with X{9) = X{0 + 1), /i(e' + 1) =/i(6') forOeTR and 



Then there exists a unique, right maximal, regular solution (/,A,/i) o/(3.3)-(3.7) 

° ° o 

with {f,X,fi)(to) = {f,X,fj,) on a time interval [to,T[ with T (E]to,oo]. If 



sup^^^i{t,e)\eelR, te[to,T[| 



< oo 



then T = 00. 



This is the content of |16, Thms. 3.1 and 6.2]. In fact this was only stated in 
the case to = 1 in that reference, but that was an arbitrary choice which makes 
no essential difference. For a regular solution all derivatives which appear in the 
system exist and are continuous by definition, cf . [|6| . 

We now establish a series of estimates which will result in an upper bound 
on /X and will therefore prove that T = oo. Similar estimates were used in for 
the Einstein- Vlasov system with Gowdy symmetry. In what follows constants 
denoted by C will be positive, may depend on the initial data and may change 
their value from line to line^ 

Firstly, integration of ( |3.5|) with respect to t and the fact that p is non- 
negative imply that 



toie- 



-, -1 



+ £) Stt 

t ^ t 



s'^p{s,9)ds 



> ■ 



t 



C~et 



, te[to,T[. 



Next we claim that 



e''+^p{t,e)de<Ct-^+', te[to,T[. 



(3.12) 
(3.13) 



A lengthy computation shows that 



dt 



p{t,d)d0-- 



2p + q-^il + ee'n 



de. 



Now q>{) and p< p so that for e = 0, 



I j\^+^p{t,e)de<~-j\^+^p{t,e)dd 
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and integrating this with respect to t yields (3.13) for e = 0. For e = — 1 we have 



2/9 + 9- 



p+p 



cW-- 



C + t 



< 



1 



c+tj^ 



e''+^{2p + q)de< 



C + tJo 



pde, 



and integrating this inequahty with respect to t yields ( |3.13 ) for e = — 1. Using 
( p3| ) and (l^ej) we find 









/ p{t,(j)da 




I 
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no /• 1 /• 1 

pie{t,T)dTda < / \^ig{t,T)\dTda 
Jo Jo 



< Ant 



oM+A 



\j{t,T)\dT<ATTt I e^'^^p{t,T)dT 



Next we show that 



< Cf, te[to,Tl de[0,l]. 



Mt,o)-x{ue)^fj^i+e^ te[t„,Tl 61 e [0,1]. 



(3.14) 
(3.15) 



To see this observe that by (gj), (|3^) and ( |3.12| ) 



d 



e"-^ ['iTTte^'^ip-p) + (ep> + l)/t] < e^-^{te^^' + l)/t 



< 



t C~et 



and integrating this inequality with respect to t yields (3.15). 

We now estimate the average of p over the interval [0, 1] which in combination 
with (3.14) will yield the desired upper bound on p: 



p{t,e)de = / p{to,e)de+ 



Pt{s,6)d9ds 



to Jo 



< C- 



f [e^''{87Ts^p + e) + l]d9ds 

J to 2s Jo 



C+-Mt/to) + 4Tr / s / e^-^e^+^pdOds + e 



to JO 



to 2s JO 



e'^'^deds 



<c+-Ht/to) + c I s 



^+'s'-'ds- 



ds 



to 



C- 



where we used ( p^ , ( |3J5| ) and (|1|). With (|1|) this imphes 

. . r e ,te[io,r[, 0e[O,l], 6 = 0, 
' ^- 1 1+lnt ,te[to,r[, ee%i], e = -i, 



(3.16) 



which by Proposition 3.2 implies T — oo. Thus we have proven: 
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Theorem 3.3 For initial data as in Proposition S.i the corresponding solution 



exists for allt€ [to,oo[ where t denotes the area radius of the surfaces of symme- 
try of the induced spacetime. The solution satisfies the estimates ( p.l3| ), (3.15) 
and (3.16). 



4 The contracting direction 

In this section we consider the Einstein- Vlasov system in the contracti ng d irec- 
tion and use conformal coordinates in which the metric takes the form (2.2). As 
in the case of areal coordinates used in the expanding direction / depends only 
on 

t, e, p\ {p' f + sin\i^x{p'^ f. 



We will write 

where 

As above, 



f = f{t,9,w,L) 
w — e^pK L-- R^{{py+smh'^x{p^f) 



p° = e-Vl + w2 + L/i?2=: e-'^ip), 



and L is a conserved quantity along particle orbits. 

In the variables introduced above the Einstein- Vlasov system can be written 
in the following form: 



dtf + ^def- 
(p) 



-■qtw-rje (p) + 



w 



L 



RHp) 



d^f^O. (4.1) 



where 



-ritRt + riBRe + — 
IK 



Rt-Re-e 



2rj 



Rte -rjtRe- Ve Rt = 47ri?e^'', 



Jim 

Pit, 9) 
qit,9) 



Rtt — Ree 



m - rjee - 



R 
1 

i?2 



Rt — Rg — e 



2') 



Rt —Re — e 



2') 



= inRe^'^p, (4.2) 
(4.3) 

inRe^'^ip-p), (4.4) 
Ane^'^ip-p-q), (4.5) 



TT 



l*00 /"OO 

P2 , , (p) fit,9,w,L)dLdw = e-^'^To„it,9), 

(^oo poo 

wfit,9,w,L)dLdw - 

-ooJO 

/•oo /'Oo 2 



TT 



e-^^n,it,9), 
R^J Jo (^^(^'^'"''^)^^'^^ = ^""'^"(^'^)' 



TT 



KP) R 



(4.6) 
(4.7) 
(4.8) 
(4.9) 
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The equations (4.2), (4.3) are the constraints and (4.4), (4.5) are the evolution 
equations; from the relation of their right hand sides with the energy momentum 
tensor it is obvious how this form of the field equations is obtained. 

Let a smooth solution of the system (4.1)-(4.9) on some time interval ]t-,to] 
be given. We want to show that if this interval is bounded and if R is bounded 
away from zero on this interval then /, i?, r] and all their derivatives are bounded 
as well, with bounds depending on the data at t = tQ and the lower bound on R. 
To this end we define auxiliary variables 



so that 



1 



The analysis which follows is modeled on the one in H, cf. also J^. 

Step 1: -hounds on R. A short computation using ([4.2|) and (4.3) shows that 



v2 



(V277C- 



1 



RrjRr 



V2R 



2V2i? 



-2%/27ri?e2''(p-j) 



observe that by ( |4.6| ) and ( [4.7D , |j|<p. Assume that R^{t,6) = for some <G 
]i_,to] and 9^TR. Then by the periodicity of R with respect to 6, 



= 0, 



= R^{t,9 + 1) < R^{t,e)eyi-p (^J^ (^V2r]^ + -^Rr'j d0' 
a contradiction. Thus Re^O on ]t^,to] x S^. Similarly, 

deRr > - (\/2?7r + ^l^^c) 
which yields the same assertion for R^- This implies that the quantity 

g°'^^d^.Rd^pR = e~^'^{-Rt^+Re^) = -]^e-^'^ R^R^ 

does not change sign. Since R is periodic in 6^ there must exist points where 
Rg = 0, hence the quantity above is negative everywhere, and by our choice of 
time direction, 

Rt>0, \Re\<Rt on ]t., to] xS\ (4.10) 



By (^), 



-Rct 



2R 



R^Rr ' 



=2») 



R 



A7rR{p-p)]>-—R^Rr; 
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observe that by (gj) and (gj) p<p. Now we fix some {t,9) G]<_,io[xIR. Then 
for se[t,to], 

^R^{s,e + t-s)^V2R^r>--^[inR{s,e + t-s)]R^{s,e + t-s). 
Integrating this differential inequality yields 

Rdt,o)<^^^^^^±^Rdto,e+t-t^ 



Similarly, 



R(t,e) 

and both estimates together imply that Rt is bou nded f rom above on ]t_,io] x S'^ 
with a bound of the desired sort. Together with ( 4.10 ) this provides bounds for 
Rt and Rg. Note that this argument shows that the spacetime gradient of R^ 
is bounded even without the assumption that R is bounded away from zero. 
Step 2: -bounds on rj. From (4.4) and (4.5) we find 



'7«r = -^i?5.-27re%. 

We fix some (i,6l) e]t_,io[xIR. Then for se[t,to]i 

-^r].{s,e + t-s)^-^-^R(is,e + t-s)~2V2'Ke^'^q{s,e + t-s). 
as Rds 

Integrating this and integrating the term containing R^ by parts yields 

R^to^e + t-to) , Ri{t,d) 



R(to,0 + t-tQ) R{t,l 



s,9 + t— s)ds. 



(4.11) 



We already know that 



/■to I 

/ R^r{s,0 + t-s)ds = -^{Rdto,0 
Jt V 2 



+ t-to)~Rdt,0)) 



is bounded, and on the other hand by (4^) the left hand side can be written as 



(s,0 + t-s)ds. 



Here the first term is bounded by Step 1, and the second and third terms are 



non-negative. Thus by (Uj), (4.6), (WJ) 



^ ° e^\{s,9 + t-s)ds< ^e^'^ip-p)is,e + t-s)ds 
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is bounded as well; recall that R is bounded away from zero by assumption. 
Thus (4.11) implies that 77^ is bounded on ]i-,io] x S^. Analogously to (4.11) 
we have 



R{to,e-t + to) R{t,e) 
+ V2^*" (2Tre^'^q-^^^ {s,e-t + s)ds (4.12) 

from which we can conclude that j]t is bounded on ]i_,to] x as well. Thus 
also rit,'i]g and rj itself are bounded there. 

Step 3: Bounds on matter quantities. We have for any tG]t_,to], 0GlR, w€ 
IR, L>0, 

f{t,0,w,L)^f{to,e{tQ,t,0,w,L),W{to,t,9,w,L),L) 
where Q{-,t,6,w,L),W{-,t,6,w,L) is the solution of the characteristic system 

,2 \ 



w = -ritw -mi {p) + TT ) + e ^"^Re 



{p)J R^P) 

of the Vlasov equation with Q{t,t,9,w,L)=9, W{t,t,9,w,L)=w. This repre- 
sentation of / implies immediately that / remains non-negative and bounded 
by its maximum at t = to- By Steps 1 and 2 the right hand side of the second 
equation in the characteristic system is linearly bounded in w, and hence if the 
w-support of / is compact initially, 

snp{\w\\{9,w,L) g supp/(t), t e]t_,to]} < 00. 

This immediately implies that p,p,j,q are bounded on ]t_,io] x S^. 
Step 4- Bounds on second order derivatives of R and rj. Steps 1, 2 and 3 
together with (O), (4.3) and (14) imply than Ree,Rte and Rtt are bounded as 
claimed on ]t_,to] x S^. The bounds on the second order derivatives of rj are a 
bit less trivial: We add the equations (4.11) and (4.12) to obtain a formula for 
r]g. When this formula is differentiated with respect to 9 there results a number 
of terms which are bounded by the previous steps and the terms 



2tt e^''— / —def(s,9 + t-s,w,L)dLdwds, 
Jt R J~ooJq kp) 

rU, _ poo poo T 

2tt / e^''— / / —def(s,9-t + s,w,L)dLdwds. 

Jt R'^J^ooJg kp) 



(4.13) 
(4.14) 



To deal with the first term we introduce the differential operators 

W = V2dr = dt-de, S^dt + ^de 

{p) 
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so that 

By the Vlasov equation 



(P) 



{?)- 



-{S-W). 



dwf- 



When this is substituted into ( 4.13| ) the resuhing term can be integrated by 
parts with respect to w, and all the terms which then appear are bounded by 
the previous steps. As to the M^-contribution, 

iWf)is,0 + t-s,w,L) = -^[f{s,9 + t-s,w,L)] 



so that the corresponding term in (i.l3) can be integrated by p arts w ith respect 



to s which again results in bounded terms. In order to deal with ( 4.14 ) we replace 
w by — w in the integral and redefine 

KP) 

the rest of the argument should then be obvious, and "qeg is seen to be bounded 
on ]i_,io] X '5'^. This way of dealing with derivatives in connection with the 
Vlasov equation was introduced for the Vlasov-Maxwell system in B. Now 



that 'qgg is bounded the same is true for r]tt by (l.f:), and -qte can be dealt with 



and (4.12) gives the necessary formula 



like rjgg: Taking the difference of (4.11) 
for 774. 

Step 5: Higher order derivatives. Via the characteristic system the C^-bounds 
on R and ry give bounds on the first order derivatives of Q{-,t,9,vu,L) and 
W{-,t,9,w,L) with respect to 9,w,L. This yields corresponding C^-bounds first 
on / and then as in Step 3 on p,p,i,q. These in turn imply C'^-bounds on R. 
The third order derivatives of 77 then have to be dealt with by repeating the ar- 
gument of Step 4, and although the details would be tedious it should be clear 
that this process can be iterated to bound any desired derivative on ]t-,io] x 
in terms of the data at t — t^ and the positive lower bound on R. 

Later we will require a slight generalization of these results. The essence 
of the first part of Step 1 above is to show that the gradient of R is timelike. 
The argument is carried out on a region which is covered by Cauchy surfaces 
of constant conformal time. However this fact is also true for any solution of 
the Einstein- Vlasov system with hyperbolic symmetry and a compact Cauchy 
surface. This follows from |19, Lemma 2.5]. Once this has been established 
the estimates in the later steps hold for any subset Z of the half-plane t<to 
provided Z is a future set. By definition this means that any future directed 
causal curve in the region t<tQ starting at a point of Z remains in Z. (For 
information on concepts such as this concerning causal structures see e.g. Q.) 
Thus if R is bounded away from zero on Z and t is bounded on Z then all the 
unknowns and their derivatives can be controlled on Z. 
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Now consider a special choice of the subset Z, namely that which is defined 
by the inequalities t\<t<tQ and 0i + io — * < ^ < ^2 — + * for some numbers B\ , 
di and t\ satisfying the inequalities B\ < 62 and ti > to ^ (l/2)(^2 — ^'i)- Suppose 
a solution of the Einstein- Vlasov system in conformal coordinates defined on 
this region is such that R is bounded away from zero. Then the functions 
defining the solution extend smoothly to the boundary oi Z at t — ti. They define 
smooth Cauchy data for the Einstein- Vlasov system. Applying the standard 
local existence theorem (without symmetry) allows the solution to be extended 
through that boundary. Repeating the construction of conformal coordinates 
in Section ^ then shows that we get an extension of the solution written in 
conformal coordinates through that boundary. 

5 Existence of an areal time coordinate 

In Section |2| it was shown that the maximal Cauchy development of initial data 
for the Einstein- Vlasov system with hyperbolic symmetry admits conformal 
coordinates in a neighbourhood of the initial hypersurface. The aim of this 
section is to show that a spacetime of this kind always admits a symmetric 
compact Cauchy surface of constant areal time. The basic strategy follows that 
of m. The fundamental object of study is a spacetime in conformal coordinates 
which develops (to the past) from initial data on the hypersurface t = 0. The 
solution will exist on a subset of the region t < of JR^. We will consider solutions 
defined on open subsets of the half-plane t<0 containing t = which are future 
sets. The union of all regions admitting solutions of this type is an open future 
set. Denote it by V. By definition i? is a positive function defined everywhere 
on V. 

Consider first the case that the past boundary of T> is empty, i.e., that 
'D—] — (X),0]xS^. Let Sp be the level set defined by the equation R{t,9)=p 
for some p less than the minimum of i? on t — Q. Since the gradient of R is 
everywhere timelike the level sets Sp are smooth spacelike submanifolds. These 
level sets can be represented in the form t = f{9), and since the coordinates 
are conformal and the level sets spacelike it follows that |/'|<1. Hence the 
Sp are compact and their pull-backs define compact spacelike hypersurfaces in 
a Cauchy development of the initial data of interest. By jsj these are Cauchy 
surfaces. As a consequence the maximal Cauchy development of those data 
contains symmetric compact Cauchy surfaces of constant areal time. 

It remains to consider the case where P has a non-empty past boundary. 
Call this boundary B. It is achronal. Hence for any given 9 there is only one t 
such that {t,9) £B. If p is a point of B then all points of B close to p are outside 
the light cone of p. Hence S is a Lipschitz curve with Lipschitz constant one. 
In other words it can be represented in the form t ~ h{9) for a function h with 
Lipschitz constant one. This is a special case of a general property of achronal 
boundaries, cf. ^, p. 187]. As a consequence B is compact. Since the gradient 
of R^ is bounded on the whole region where the solution is defined (as shown 
in Section 0) it follows that R^ is uniformly continuous. Hence the function R^ 
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extends to a continuous function on the closure 2? of 2?. Its square root i? is a 
continuous extension of i? to P. It will now be proved that -R = on B. From 
this point on it is possible to argue as above to obtain the existence of a Cauchy 
surface on which R is constant. 

Consider first a point 9o where h has a local maximum and let to = h{6o)- 
Then there exists e > such that for any t greater than to the points {t,9) with 
00 — e <9 < Oo + e lie in V. If R{tQ,9o) > then the solution has a smooth limit 
at the points {to,0) with 9o — e<9 <0o + e provided e is small enough. This 
gives Cauchy data at t = to- The corresponding local solution of the evolution 
equations can be used to extend the original solution, thus contradicting the 
assumed maximality of T). It follows that in fact R vanishes at any point of B 
where h has a local maximum. 

Let Vi be the intersection of V with the hypcrsurface where the time coor- 
dinate takes the value t. It is an open subset of S^. It is either the whole of 
S^, the empty set or a disjoint union of open intervals which are its connected 
components. Consider now a point (fo,^o) of B which is an endpoint of a com- 
ponent of T>tQ. Suppose further that R{to,9o) >0. Without loss of generality 
we can assume that is a right endpoint of a component; the argument for a 
left endpoint is strictly analogous. For some e > the solution extends smoothly 
to the points {to, 9) with 9o — e<9<9o- This can be extended to smooth data 
at t = to on the interval 9o — €<9<9o + e- There is a corresponding local solu- 
tion. Using the domain of dependence we see that this new solution extends the 
original solution at those points outside the past light cone of the point {to,9o)- 
In this way we obtain an extension of the original solution and a contradiction 
unless the part of B immediately to the left of 6*0 is a null curve. If Bi is the set 
of points of B which are endpoints of components of some Vt and which have 
.R>0 then it has now been shown that Bi is a union of null segments. Return- 
ing to the point (to,9o) we can ask what happens to the interval of which it is 
an endpoint as t increases. The endpoint could vanish if the interval coalesces 
with another while R remains positive. However that would mean that 9o would 
have to be a local maximum of h, a case already excluded. Hence as we move 
to the right in 9 the corresponding point of B must continue to be an endpoint. 
Eventually R must become zero. But in that case the future-pointing character 
of the gradient of R can be used to obtain a contradiction. It can be concluded 
that Bi is empty. 

It remains to consider the case that there is a point {to,9o) of B which is 
not an endpoint of a component of some P*. In this case {t,9o) belongs to T> 
for every t>to- Let {a{t),b{t)) be the component of Vt containing ^o- The 
functions a{t) and b{t) are decreasing and increasing, respectively. Let ao and 
bo be their limits as t— >to from above. If ao < bo then {to,9o) corresponds to a 
local maximum of / and so ^ = there. Otherwise it is a limit of the points 
{t,b{t)) which are right endpoints of components of Vt, and so by continuity 
.R = 0. 

A natural question arising in this context is whether a spacetime which exists 
globally in one time direction with respect to a conformal time coordinate cannot 
be extended in that direction to a strictly larger globally hyperbolic spacetime. 
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It will be shown that this is impossible. As an aid in doing this we consider the 
action of the symmetry group on M. We know exactly what the action looks 
like near the initial hypersurface. There is a two-dimensional subgroup H oi G 
which acts transitively on the orbits of the action of G. Corresponding to H 
there are two Killing vector fields which are spacelike and linearly independent 
near the initial hypersurface. We claim that they are spacelike and linearly 
independent on the whole of M, so that the orbits of the action of H are 
everywhere spacelike and two-dimensional. To show this it suffices to show that 
a Killing vector generated by the action of H can neither vanish nor become 
null anywhere on M, cf. |0, Prop. 2.4]. Suppose that a Killing vector of this 
kind vanishes at a point p. Without loss of generality we can assume that p lies 
to the future of the initial hypersurface and that the Killing vector is spacelike 
and non-zero on the chronological past of p. The integral curve of the Killing 
vector through p is null and if we extend it maximally it must either reach the 
initial hypersurface or tend to a zero of the Killing vector. The first of these 
possibilities contradicts the fact that the Killing vector is spacelike on the initial 
hypersurface. In the second case we see that the Killing vector has a zero to 
the future of the initial surface. 

Now it will be shown that the existence of zero of the Killing vector at a 
point p to the future of the initial hypersurface leads to a contradiction. The 
flow generated by the Killing vector leaves p invariant. Its linearization maps 
the space of null vectors at p to itself and must have a fixed point. In other words 
the linearization leaves a null direction invariant. The null geodesic with this 
initial direction is then also invariant. It must intersect the initial hypersurface, 
and the point of intersection is invariant under the flow. This gives the desired 
contradiction. 

Now we come back to the question of inextendibility. Suppose that the 
spacetime exists for all conformal time in one time direction, without loss of 
generality the future. Suppose that there is a point p belonging to a future 
development of the initial data but not in the region covered by the conformal 
coordinates. It lies on an orbit of H. Let 71 and 72 be the two null geodesies 
through p orthogonal to the orbit. They are orthogonal to all orbits they meet. 
They must enter the region covered by conformal coordinates and then they 
coincide with the curves with t = ±9 + C for some constant C and constant 
values of the other coordinates. Without loss of generality we can assume that 
p is on the boundary of the region covered by the conformal coordinates. In the 
approach to p the geodesic 72 intersects the geodesic 71 infinitely many times. 
This contradicts the strong causality and hence the global hyperbolicity of the 
extension. The conclusion is that an interval of conformal time which is infinite 
in one time direction proves that the maximal Cauchy development is exhausted 
in that time direction. 

There is also another criterion for inextendibility which can conveniently be 
discussed here. It is related to an argument given in j^]. Consider the action 
of the group H on M introduced above and the corresponding Killing vectors. 
Since they are always linearly independent with spacelike orbits, the 2x2 matrix 
of inner products of these vectors always has positive determinant A. Suppose 
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that a solution is given which exists globally towards the future in areal time. 
It will now be shown that it has no proper globally hyperbolic extension to the 
future. If there was such an extension there would be a point p of the extension 
not belonging to the original spacetime. Let 71 be a null geodesic through p as 
above. We may assume without loss of generality that it immediately enters the 
original spacetime on leaving p towards the past. In the region covered by the 
areal time coordinate it has constant values of the coordinates x and y. Along 
the part of 71 belonging to the original spacetime A is proportional to i?^. As 
a point on 71 approaches p the areal time coordinate t tends to infinity, which 
means that A tends to infinity. But A is a smooth function at the point p. 
Hence no point p of this kind can exist. 



6 Main results 

In this section the analytical and geometrical information obtained in previous 
sections is combined to obtain the main results of the paper. 

Theorem 6.1 Let {AI,gaf}, f) be the maximal globally hyperbolic development 
of initial data for the Einstein- Vlasov system with hyperbolic symmetry. Then 
M can be covered by symmetric compact hypersurfaces of constant area radius. 
The area radius of these hypersurfaces takes all values in the range ]i?o,oo[ where 
Rq is a non-negative number. 

The quantity Rq depends on the solution, but there is no good understanding 
of how many solutions have i?o 7^ 0. 

Proof of Theorem \6.^ . For a spacetime satisfying the hypotheses of the theorem 
we know from Section ^ that a conformal coordinate system can be introduced 
on a neighbourhood of the initial hypersurface 5*0 corresponding to the original 
data. By the results of Section || it follows that this region can be extended to 
the past so as to include a Cauchy surface Sa of constant areal time. Moreover, 
either the conformal time coordinate extends to all negative values, or R tends 
to zero as the past boundary B of the region covered by conformal coordinates 
is approached. In the first of these cases the region covered by the conformal 
time coordinate includes the entire past of the initial hypersurface in the max- 
imal Cauchy development, as shown in Section ||. Also the past of Sa in that 
region admits a foliation by hypersurfaces of constant R. In that region we can 
transform to areal coordinates. For we can choose the spatial coordinate 6 so 
that its coordinate lines in Q are orthogonal to that foliation. In the second 
case (where R tends to zero on a boundary B) the past of Sa is also covered 
by areal coordinates. It exhausts the past of Sa in the maximal Cauchy devel- 
opment, as will now be shown. A well-known argument related to Hawking's 
singularity theorem shows that it is enough to check that the mean curvature of 
the foliation tends uniformly to infinity as R tends to zero. The mean curvature 
is given in areal coordinates by |l^ 

trfc = -e-''(A + 2/t) (6.1) 
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Substituting the field equation for A and using the fact that p > gives the 
inequality 

trfc<-(l/2t)(3e-^ + e^)<-\/3/t (6.2) 

and hence the desired result. Thus in both cases the entire past of Sq in the 
maximal Cauchy development is covered. It follows from the results of Section |^ 
that the spacetime and the areal time coordinate can be extended so that the 
time coordinate covers the interval ]i?o,oo[. It can then be concluded by the 
argument at the end of Section |^ that the entire future of 5*0 in the maximal 
Cauchy development is covered. 

Theorem 6.2 Let {M,gap, f) be the maximal globally hyperbolic development 
of initial data for the Einstein- Vlasov system with hyperbolic symmetry. Then 
M can be covered by compact hypersurfaces of constant mean curvature. The 
mean curvature of these hypersurfaces takes all values in the range ] — oo,0[. 



Proof. It follows from ( p. 2D that there are Cauchy surfaces with everywhere 
negative mean curvature. Under these circumstances it was shown by Henkel 
[ pi] that the initial singularity is a crushing singularity and thus a neighbour- 
hood of it can be foliated by CMC hypersurfaces. Next the statement about the 
range of the CMC time coordinate will be proved. The important observation 
is that the argument used to prove the corresponding statement in [plf extends 
straightforwardly to the case of hyperbolic symmetry. It is more powerful than 
the arguments used to extend the CMC foliation in jl9|. It remains to see 
that the CMC foliation covers the entire future of the initial hypersurface. It 
is enough to show that if p is any point of the spacetime there is a compact 
CMC hypersurface which contains p in its past. Let Si be the Cauchy surface 



of constant areal time passing through p. Equation (6.2) shows that the mean 
curvature of 5*1 is strictly negative. Hence it has a maximum value Hi < 0. Let 
5*2 be the compact CMC hypersurface with mean curvature Hi/2. Then the 
infimum of the mean curvature of ^2 is greater than the supremum of the mean 
curvature of 5*1 and a standard argument JTsI shows that ^2 is strictly to the 
future of Si. Hence p is in the past of 5*2, as required. 

Remark. In the case of Gowdy symmetry the mean curvature of the hypersur- 
faces of constant areal time is negative and the corresponding argument applies, 
thus showing that the CMC foliation covers the spacetime in that case too. 

The considerations of this paper were confined to the case where the matter 
content of spacetime is described by the Vlasov equation, but in fact many of the 
arguments do not depend on the details of the matter model. It would be worth 
to investigate which parts of the conclusions extend to which matter models. 
Some relevant estimates for scalar fields, wave maps and electromagnetic fields 
can be found in @, |2l| and Q. 



7 Possible further developments 

In this paper we have seen how it is possible to get rather complete information 
on the existence of global geometrically defined time coordinates in spacetimes 
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with hyperbolic symmetry and, more generally, with surface symmetry. The 
motivation for being interested in this question was that its answer could provide 
a tool for understanding the global geometry of these spacetimes. To do so it 
is necessary to understand something about the asymptotic behaviour of the 
solutions as expressed with respect to these time coordinates. The kind of 
questions which we would like to answer are the following. Is the spacetime 
geodesically complete towards the future? Is there a curvature singularity in 
the past? Is the singularity in the past velocity dominated? Does the solution 
become homogeneous in the future in some sense? Can we obtain detailed 
asymptotic expansions for the solution in these regimes? 

The one strong and rather general result concerning questions like this which 
is available is the following theorem of |0 . If there is a foliation by Cauchy sur- 
faces of constant R such that R approaches zero as the singularity is approached 
and if, in the hyperbolic case, the Hawking mass is positive near the singularity 
then the Kretschmann scalar blows up as the singularity is approached. That 
some restriction is required is shown by the pseudo-Schwarzschild solution (cf. 
[|9)). For negative values of the mass parameter this solution, which has hy- 
perbolic symmetry, has negative Hawking mass, and R does not approach zero 
at the boundary of the maximal Cauchy development. The Kretschmann scalar 
remains finite near that boundary. A guess, which is not contradicted by any 
known results, is that the pseudo-Schwarzschild solutions are the only solutions 
of the Einstein- Vlasov system with hyperbolic symmetry which do not have a 
singularity where the Kretschmann scalar blows up. It would be very interesting 
to have a rigidity theorem of this kind. 

Apparently the only other source of information about the asymptotics of 
surface symmetric solutions of the Einstein- Vlasov system comes from the study 
of spatially homogeneous solutions. The spatially homogeneous solutions with 
spherical, plane and hyperbolic symmetry are the spacetimes of Kantowski- 
Sachs, Bianchi I and Bianchi III types respectively. There are exceptional cases 
(all vacuum) which are the Schwarzschild solution (appropriately identified), 
the flat Kasner solution and the pseudo-Schwarzschild solution. These were 
presented in Appendix B to ||l^. All other solutions are such that curvature 
invariants blow up in both time directions (Kantowski-Sachs) or that curva- 
ture invariants blow up in one time direction and the spacetime is geodesically 
complete in the other time direction (Bianchi types I and HI) po| . 

More detailed results on the asymptotics have been obtained for Bianchi 
types I and III |Q . Translating the theorems obtained there into the notation 
of this paper leads to the following results. All non- vacuum solutions of either of 
these two Bianchi types are such that i? ^ as the past boundary of the maximal 
Cauchy development is approached. Generically /i~ (l/2)lnt and A~ — (l/2)lnt 
as where we only indicate the leading term in an asymptotic expansion. 

There is a smaller class of solutions (not explicitly known) where /i^(l/2)lnt 
and A ~ const, as >0. Finally there is an even smaller class where /i~lnt and 
A~lnt as t^O. In the expanding direction all solutions of both Bianchi types 
are geodesically complete. For type I we have /i ^ (1/2) \nt and A ~ \nt as t ^ oo. 
The behaviour of the type III solutions was not completely determined in p2|. 
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It was shown that (in the notation of the present paper) A is increasing for large 
t but probably it grows slower than any positive multiple of Int. 

Clearly the next step is to extend these results on asymptotic behaviour 
to the inhomogeneous case. It would be convenient if the inhomogeneous solu- 
tions behaved essentially like the homogeneous ones in their asymptotic regimes. 
There are, however, phenomena which may prevent this. These are the possi- 
bility of the occurrence of analogues of the spikes found in Gowdy spacetimes 
[ p3| in the contracting direction and the Jeans instability Q in the expanding 
direction. 
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